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Abstract 

When  a  high  energy  solid  propellant  or 
explosive  becomes  porous,  either  by  design  or  by 
accidental  damage,  it  becomes  more  sensitive  to 
shock-to-detonation  transition  (5DT).  Stress 
waves  propagating  ahead  of  the  convective  flame 
front  cause  pore  collapse  which  creates  a 
confined  bulk  reaction  zone.  With  sufficient 
pressur 1 zat ion  In  the  reaction  zone,  the  stress 
waves  coalesce  to  form  a  shock.  This  paper 
describes  an  attempt  to  model  the  shock  formation 
In  porous  solids  by  coupling  the  gas 
pressurization  to  the  solid  mechanics  of  pore 
collapse  and  shock  formation.  Concepts  are 
presented  which  Indicate  why  porous  materials  are 
more  sensitive  to  shock  Initiation  than  solids; 
mainly  It  Is  due  to  a  greater  bulk  energy  and 
formation  of  localized  hot  spots. 

1.  Introduction 


granulated  compressed  MMX  explosive,  that 
detonation  may  be  a  result  cf  shock  Initiation. 
To  explain  Campbell’s  results,  a  new  mechanism 
describing  the  detonation  of  porous  propellants 
was  given.  Compression  waves  originating  at  the 
burning  front  propagate  ahead  causing  compaction 
of  the  unreacted  porous  mixture  resulting  In  the 
formation  of  a  solid  plug  thereby  eliminating 
further  convective  flame  propagation.  The 
compression  waves  continually  generated  at  the 
burning  front  propagate  through  the  solid  plug. 
Coalescence  of  these  compression  waves  form  a 
shock  wave  which  Initiates  the  explosive 
resulting  in  a  detonation. 

Work  done  by  Macek  [4]  and  Tarver  et  al.[5] 
used  the  method  of  characterl st Ics  to  predict 
shock  development  in  homogeneous  solid  explosives 
from  confined  end  burning.  In  both  studies  a 
modified  Talt  equation  was  used  as  the  solid 
material  state  relation. 


Advanced  formulations  for  solid  propellants 
continue  to  be  more  energetic  and  the  possibility 
o  deflagration  to  detonation  transition  (DDT)  is 
real,  especially  If  the  propellant  motor  grain 
becomes  porous  due  to  damage  by  accidental 
stimuli.  In  a  recent  transient  analysis  [1]  for 
DDT  In  a  porous  propellant  bed.  It  was  shown  that 
steady  detonation  characterl sties  (CJ  pressure 
CJ  temperature,  detonation  velocity)  were 
predicted  without  requiring  the  precursor  shock 
common  to  most  detonations  In  homogeneous 
materials.  This  was  referred  to  by  some  as  the 
convective  burn"  model.  In  Ref.  [1]  It  was 
shovffi  that  the  high  $urface-to- volume  ratio  of 
very  small  (  l/10irtn  diameter)  propellant 
partldes,  high  burning  rate-pressure  Indices 
(0.9  to  1.0),  and  confinement  of  product  gases 
(due  to  the  tight  packing  and  enclosing  walls) 
combined  to  cause  a  rapid  enough  gas  generation 
rate  to  build  and  sustain  pressures  equivalent  to 
those  observed  In  granular  bed  experiments 
L?,3J.  However,  as  stated  above,  the  model 
analysed  In  Ref,  [1]  did  not  attempt  to  treat  the 
stress  waves  propagating  Into  the  unbumed  porous 
material.  It  Is  the  purpose  of  this  paper  to 
analyse  the  propagation  and  shock  formation  of 
these  compression  waves  generated  at  the 
products-reactants  Interface. 


r  ..  ?iarI?  r  *  receflt  set  of  exper Indents ,  A.W. 
Campbell  13J  has  sho*wi  by  using  a  similar  test 
arrangement  to  that  useo  In  [2j,  but  with 


Macek  used  the  method  of  characterl sties  to 
theoretically  predict  his  observed  experimental 
results  [41.  By  assuming  an  exponential  type 
pressure  rise  1r  the  product  gases,  Macek  tr»ced 
right  running  characterl st 1 cs  generated  by  * 
pressure  input  and  graphically  showed  the  ,r .>* 
to  coalesce  at  approximately  the  distant- 
observed  In  his  experiments.  The  coalesce  :e  ,f 
the  characteristics  Is  Interpreted  as  the 
formation  of  a  shock  wave  which  Is  assume 
cause  initiation  of  the  explosive  resuitln  *n  a 
detonation  propagating  forward  and  i  relation 
propagating  backwards. 


The  work  presented  in  this  pap*r  will  show 
that  the  Talt-equatlon  assumption  i\ade  in  [4]  and 
[5]  was  Incorrect  for  this  particular  uompressl on 
process  and  a  more  general  Mie-Grunelsen  equation 
of  stite  (E.O.S. )  should  be  used. 

In  fact,  a  method  of  character  1 s t tes 
calculation  Is  Impossible  with  the  Mie-Grunelsen 
state  equation,  so  the  system  of  conservation 
equations  were  solved  numerically  by  a  Lagrangtan 
finite  difference  code.  Homogeneous  solids  with 
known  properties  were  first  computed  to  test  the 
validity  of  the  code,  followed  by  calculations 
for  heterogeneous  or  porous  materials.  The 
purpose  of  this  work  Is  to  eventually  incorporate 
reactive  chemistry  and  hot  spot  mechanisms  into 
the  model  and  go  beyond  the  shock  formation  to 
predict  DOT. 


r?op£ll!ntcritS:  Sh0Ck  W4y*S  40d  Deto,,ators;  Active  Flows;  Properties  of 
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Shown  in  Figure  1  is  a  schematic 
representation  of  what  we  envision  as  the  six 
part  DSDT  process.  Superimposed  on  each  section 
of  the  figure  Is  a  line  representing  the  qas 
porosity  (gas  volume/total  volume)  as  a  function 
of  bed  location.  A  value  of  +  equal  to  unity 
represents  a  zone  of  all  gas  while  *  equal  ta 
zero  indicates  a  homogeneous  solid. 


Portion  A  represents  the  initial  blast 
caused  by  the  ignitor  used  to  Initiate  OSOT  in 
the  bed.  In  the  second  portion*  the  convectlvely 
burning  porous  bed  Is  shown  to  be  pressurizing 
the  region  behind  the  ignition  front.  This  Is  a 
result  of  the  product  gases  being  confined  by  the 
walls.  The  porosity  variation  for  this  time 
Indicates  an  all  gas  makeup  in  this  particular 
zone.  Actually,  there  is  a  finite  reaction  zone 
where  the  propellant  is  decomposing  into  the 
product  gases.  Reference  [1]  showed  that  if  the 
reactive  material  consisted  of  rm  size  particles 
this  reaction  zone  may  be  very  large.  This 
results  In  a  slower  pressurization  rate  of  the 
gas  cavity.  Illustrated  In  the  last  four  parts 
of  Figure  I  Is  what  we  refer  to  as  the  second 
possible  branch  of  OSOT.  the  shock  compression 

?;£*  POre  col1jP«.  Plug  formation 
and  $h°ck  formation  are  envisioned  concepts  which 
will  be  expanded  upon  In  the  following  sections 
of  this  paper. 


2.  Governing  Equations 

For  the  hydrodynamic  analysis  of  a 
continuous  material  with  a  moving  boundary  the 
governing  equations  are  easier  to  work  with  If 
expressed  In  the  lagranglan  or  material  form  as 
opposed  to  the  Eulerlan  form.  The  conservation 
of  mass,  momentum  and  energy  for  one-dlmenslonal 
unsteady  flow  are  written  as 


v  ux 

(1) 

**  Px 

(2) 

-P  vt 

(3) 

In  the  above  expressions  v  represents 
specific  volume,  u,  particle  velocity,  e. 
specific  Internal  energy  and  P.  the  total  stress 
acting  on  the  particle.  The  subscripts  x  and  t 
indicate  partial  derivatives  with  respecT  to  **“ 
distance  and  time,  respectl vely,  evaluated  at  a 
constant  particle  location. 


For  comparison  the  same  conservation 
equations  written  in  Eulerlan  or  spatial  form  are 

et  ♦  fcu)x  •  0  (4) 

(pu)t  ♦  feuu  ♦  P)x  -  0  (5) 

(pe)t  *  (pue-+  Pu)^  *  0  (6) 
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fig.  1  Six  part  representation  of  OSDT 

(A)  Ignition  of  porous  material,  (B) 
Convective  flame  spreading  with  pressure 
rise  In  gas  zone.  (C)  Pore  collapse, 

(D)  Plug  formation,  (E )  Shock  formation, 
(F)  Shock  Initiation  of  porous  material 


*  OSOT  •  Deflagration  to  ^hock  to  Detonation  transition 


2 


derivative. 


Here  the  partial  derivatives  are  evaluated  in  a 
fixed  reference  frame  and  the  material  particles 
are  In  motion  relative  to  that  frame,  thus, 
giving  the  x-dlrectlon  flux  terms. 

A  fourth  expression,  an  equation  of  state 
for  the  material,  forms  a  mathematical  relation 
between  pressure,  specific  volume  and  specific 
energy.  The  particular  form  of  the  state 
equation  used  will  be  expanded  on  in  the 
following  section  and  for  the  time  will  be 
expresses  .  the  general  form  as 

?  *  (7) 

Equations  (1-3)  and  (7)  represent  the  full  system 
of  equations  necessary  to  describe  the  flow 
phenomena  In  a  homogeneous  material.  More 
specifically  they  are  four  Independent  equations 
In  four  unknowns:  u,  v,  P#  and  e.  Since  thermal 
conduction  Is  neglected  In  oi  *  calculations, 
temperature  does  not  appear  In  the  governing  flow 
equations.  For  this  reason  temperature  can  be 
treated  as  an  "auxiliary"  relation  and  will  be 
dealt  with  later  in  the  paper. 

When  describing  stress  wave  propagation 
through  porous  materials,  an  additional  variable 
appears,  a  ,  the  porosity. 

a  =  v/v$  (where  a  >  1 )  (8) 

For  a  given  volume  v,  the  volume  occupied  by 
homogeneous  solid  Is  represented  by  v  .  The 
difference  between  these  two  volumes  is  treated 
as  a  vacuum  for  these  calculat Ions.  Other 
researchers  [6]  have  modeled  this  void  volume  as 
a  gas  which  Is  Iseotroplcally  compressed  as  the 
voids  collapse  and  thus  act  to  resist  the 
compaction -col lapse  process. 

With  the  addition  of  a  fifth  Independent 
variable  namely,  a  ,  an  additional  equation  Is 
necessary. 


r(v$)  =  vs(aP/de)y  (io) 

Evaluating  equation  (9)  at  ambient  conditions  and 
making  use  of  several  thermodynamic  relations 
(See  Reference  [6])  gives  the  constant. 


r'vso' u  ro  *  ~c  on 

V 

where  8  is  the  coefficient  of  thermal  expansion; 
cQ#  Is  the  ambient  sound  velocity;  c  is  the 
specific  heat  at  constant  volume.  Values  of 
these  parameters  for  the  material  HMX  are  given 
in  Table  1. 


Table  1.  HMX  Properties 


Initial  Specific  Volume 
Initial  Density 
Ambient  Sound  Velocity 
U  -  u  Slope 
Initial  Yield  Stress 
Initial  Shear  Modulus 
Specific  Heat  (constant 


vSQ  *  0.526  cc/g 

p  so  *  l-9  9/cc 
c  *  .267  cm/ysec 
S  •  2. 60 
Y  *  0.0517  GPa 
G  -  3.516  GPa 
volume)  C  -  10.  U/g°K 


A  simplifying  approxlmatl on  made  for  the 
Grunelsen  coefficient  Is 

r(vs)  p 

T  *  foV  constant.  (12) 

Equation  (9)  Is  convenient  for  shock 
development  problems,  since  It  gives  the  pressure 
at  any  thermodynamic  state  as  the  deviation  , 
aP  *  r ( v  )  (e-e^J/v  ,  from  the  volume 
dependent  reference  pressure,  PH5.  The 
subscripts  H$  In  the  above  equations  Indicate 
that  the  reference  state  chosen  Is  the  solid 
material  shock  Hugonlot. 


Carroll  &  Holt  [6]  have  done  extensive  work 
In  studying  and  mathematical ly  modeling  the 
collapse  of  porous  materials.  Our  wo^k  employs 
their  moH«i  for  static  pore  collapse,  relating 
the  stress  on  the  material  to  the  void  volume. 
For  the  time  being  it  will  be  expressed 
as,  P  •  P(  o  )•  A  more  detailed  description  of 
the  three-part  model  will  be  presented  later. 


Experimental  shock  velocity-particle 
velocity  data  indicates  a  linear  rlt  Is  valid  for 
volume  changes  encountered  In  our  analysis. 
Reference  [8]  gives  for  HMX , 

U${cm/Msec)  -  A  ♦  S  up(cm/gsec)  (13) 


where  A  -  0.267  and  S  •  2.6  . 


3.  Equations  of  State 
Homogeneous  Solids 

Most  of  the  work  presented  in  this  paper 
made  use  of  a  Mle-Grunel sen  equation  of  state  for 
homogeneous  solids  and  a  modified  form  of  It  for 
the  porous  material  calculations.  The  general 
form  of  the  Mle-Grunel sen  equation  is 

P(V*)  ■  W  +  («-«H$)r(vJ)  (9) 

V^ 

In  equation  (9),PH$  and  eHV  the  Hugonlot 
pressure  and  Hugonfot  energy,  are  used  as  the 
reference  state  and  r  (v  )  Is  the  Grunelsen 
coefficient  representing  the  thermodynamic 


In  the  limit  of  an  Inf Ini tesslmal  disturbance  In 
particle  velocity,  u  *  0,  It  becomes  obvious 
that  the  constant  A  In  equation  (13)  rtpre«ants 
the  adiabatic  sound  velocity,  c  . 

J  o* 

Combining  equation  (13)  with  the 
conservation  of  momentum  and  conservation  of 
energy  across  the  shock  discontinuity ,  an 
expression  for  the  Hugonlot  or  reference  pressure 
used  in  equation  (9)  Is  obtained 


HS 


<*,> 


»«oco*  R 

[l-SR]2 


(14) 


where  R  -  and  0  Is  the  initial  solid 

density.  so 
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Past  work  In  DOT  by  Macek  [4]  and  Tarver  et. 
a! .  [5]  used  a  modified  Taft  equation  as  their 
constitutive  relation,  namely. 


P(vs) 


p  soco 


(15) 


isentrope  AD6CA  represents  the  reversible  wor* 
done  in  an  equivalent  isentropic  compression  and 
the  remaining  slice,  A0BAt  represents  the 
Irreversible  entropy  Increase  associated  with  the 
shock  di scont i nu 1 ty .  A  qualitative  comparison 
will  he  made  when  discussing  shock  waves  In 
porous  material . 


This  particular  state  equation  was  used  with  the 
supposition  that  It  represented  isentroplc 
conpressl on  throughout  the  entire  loading 
process.  In  both  studies,  because  of  the 
simplicity  of  equation  (15)  and  the  Isentroplc 
assumption,  the  method  of  characteristics  was 
easily  applied  to  the  governing  hyperbolic 
equations  and  characterl sties  (or  distance-time 
traces  of  constant  stress)  could  be  drawn  for 
various  stress  loads  and  upstream  shock 
coalescence  predicted.  Appendix  A  includes  a 
brief  description  of  the  method  of 
characterl sties  as  applied  to  one-dimensional 
transient  flow  In  a  homogeneous  material. 

In  order  to  determine  whether  the  Talt 
equation  assumption  was  Indeed  rcoresentl t 1 ve  of 
an  Isentroplc  process  up  to  pressures  of  4GPa, 
(40  kbar),  we  plotted  In  Figure  2  the  Talt 
equation  along  with  the  shock  Hugonlot,  and  an 
Isentrope  (S  -  SQ)  for  HMX. 
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Fig,  2  Comparison  of  HMX  Shock  Huqonlot, 

Isentrope  (T/v  ■  constant^  and  Tait 
equation  of  state 

Figure  (2)  Indicates  that  the  Talt 
expression  only  models  Isentroplc  compression  for 
pressures  less  than  approximately  0,5  GPa 
(■  5  kbar  ■  72,700  pslg)  and  becomes  more 
Incorrect  the  larger  the  stress.  Also  depicted 
In  Figure  (2)  Is  the  small  difference  between  the 
shock  Hugonlot  and  Isentrope,  This  is  common  for 
solids  when  the  pressure  is  less  than 
approximately  5  GPa  [7], 

By  shock  compressing  a  solid  from  ambient 
conditions  to  the  state  labeled  B  on  Figure  (2), 
an  amount  of  Internal  energy  equal  to  the 
triangular  area  ABCA  (e«  -  eQ  •  P»  (vA  -  vR)/2) 
is  deposited  Into  the  shocked  state  (state  B). 

Of  this  total  energy,  the  area  under  the 


Porous  Materials 


The  state  equation  for  the  porous  material 
calculations  Is  the  sane  Mie-Grunei sen  equation 
used  for  homogeneous  sol  ids  with  the  introduction 
of  a  ,  the  "porosity**, 

p  (v/a.  e)  .  PHS(v/a)  *  <16> 

where  v  and  v  represent  the  spnciflc  volume  and 
the  Initial  specific  volume  of  tne  mixture. 
Equation  (16)  gives  the  pressure  for  a  porous 
material  as  a  function  of  the  energy  of  that 
material  and  a  reference  pressure,  pHS(v/  a  ). 
Since  v$  -  v/  a  ,  this  reference  curve  Is  the 
same  as  that  for  the  homogeneous  solid,  i.e,  the 
solid  material  Hugonlot, 

With  the  introduction  of  the  variable  a  , 
porosity,  an  additional  constitutive  relation  is 
needed  to  fully  describe  the  material  state.  The 
following  section  introduces  the  pressure- 
porosity  (or  P-  a  ).1aw  used. 

4,  Pore  Collapse 

A  constitutive  relation,  refered  to  as  the 
"P-  a  law,"  is  needed  when  the  conservation 
equations  are  solved  for  a  porous  solid.  The  law 
relates  an  externally  applied  pressure  to  the 
extent  of  compaction  or  pore  collapse, 
represented  by  a  *  The  form  of  the  P-  a  law 
used  was  developed  by  Carrol  and  Holt  [6],  In 
the  model,  the  porous  matrix  Is  treated  as  a 
hollow  sphere  where  the  Inner  and  outer  radii  are 
chosen  such  that  pore  size  and  overall  porosity 
of  the  porous  material  being  modeled  are 
accurately  represented. 

Pore  collapse  occurs  In  three  phases:  (1) 
elastic  phase,  where  the  solid  elasticly  deforms: 
(2)  elastlc-plast ic  phase,  where  plastic 
deformation  begins  at  the  inner  radius  and 
progresses  outward  until  plastic  deformation 
begins  at  the  outer  radius:  and  (3),  the  plastic 
phase,  where  plastic  deformation  occurs 
throughout.  The  P-  a  relations  for  the  three 
phases  of  compaction  and  the  appropriate  range 
which  each  applies  are  given  by. 


elastic  phase: 


ao>-  “  >  ai 


AG  (a  -  a) 

p  ‘  irfriT 

elastic-plastic  phase: 


i j  >.  a  > 


7  ?G(a  -a) 

P*1  *  lnf“TcT-“iyi 


(17*) 


(17b) 
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*  *  °  °r  U  ~  <c-<0>  (A-’) 

The  C  characteristics  generated  at  the  left 
end.  can  now  be  plotted  on  the  x-t  p’ane  at  my 
tlpne  for  a  given  P(t)  relation,  A  summary  of  the 

!$  9lven  b(?,ow  for  4  general  pressure 
rate  P ( t )  Input, 

Table  A- 1 

logic  To  Use  Hethod  of  Characterl st Ics 
For  Shock  Compression 


Step  f  Prccess 
"K  it  time  t,  P*P (t ) 


2. 

3. 

4. 

5. 


calculate  p/Pq  from  Tait  E.O.S.,  *q. 

calculate  the  sound  speed  from  Eq.  A.o 

calculate  the  velocity,  u,  from 
tq.  A. 9 

determine  the  position  of  the  left 
boundary  from 

V'J  u  *  • 

using  equatfon  A.9,  A.6.  A.5  and  the 
KCty  input,  u  can  be  "xpressed  In 
terms  of  t  only 

2c  h*1 

X  •/*  [— 2  ra.d(t)  IT  ,. 
o  o  ‘  n-1  a  2  dt 

p0  C0 

calculate  the  slope  of  the 
characteristic  dx/dt  ■  u  ♦  c 

plot  on  the  x-t  plane 
*(t)»  xq+  (dx/dt )t 


IK  Appendix  B 

TEHPfRAnjPE  INCREASE  DUE  TQ  SHOCK  COMPRESSION 

.  „  The  ^thod  used  for  determining  temperature 

Increases  due  to  shock  waves  Is  that  of  P 

a°d  Ra1zer  W  »nd  Involves  the 

«^ItJ°n1n9  °f  ener5ty  and  Pre5Sure  Into  three 
separate  components. 

first  component  of  pressure  Is  due  to 
th*  Interatomic  potential  and  Is  referred  to  as 
**  *l4St  c  or  te™.  At  eero  degrees  Ke?v1n 

reper'or^attract^  pretsu?'  the  ‘tews  neither 
1  °°e  another  at  this  point 
the  pressure  Is  zero.  The  volume 

occupied  by  the  solid  at  this  point  V 
slightly  less  than  the  volume  occupied^  the 
same  mass  of  material  it  room  temperature 

SToVS  tthrffla'!«/'pands  whw  bMted 

Fruw  u  *  to  JUO  K ,  the  solid  actually  a 

fo?«  V«*MblMtPreS^r?  (attr4Ct1v*  Interatomic 
caressed  !??  conditions.  As  the  material  Is 
iH  S'lul  er  ,sentrop1cally  or  Irreversibly 
I  *  process,  the  atoms  repel  one 

Sre«ure*1n  »  ?r  co"tr1hute  to  the  elastic 
pressure  in  a  positive  sense.  The  elastic 


component  cf  the  total  energy  is  related  to  the 
elastic  pressure  through  the  wort-enerqv 
relation, 


Pc  dv  =  -dec 


(B.l) 

where  the  subscript  c  is  In  reference  to  the 
elastic  or  cold  component  of  the  total  Internal 
energy,  e,  or  total  pressure.  P.  The  name  cold 
energy  or  pressure  has  significance  In  that  the 

ho  tK'Cn^eSSUr!'VOlume  reldt,on  can  be  shown  to 
be  the  0  K  isotherm  and  also  the  S  *  0  isentrope. 

The  second  and  third  components  of  the 
energy  are  due  to  the  atomic  and  electronic 
motions,  respectively.  The  elctronic  terms,  e 
and  Pe,  are  only  of  significance  above  e 

approximately  10.000'K  and  hence  outside  the 
regime  of  the  processes  encountered  here. 

»h  i/be  thermal  component  of  energy  represents 
the  harmonic  oscillations  of  the  atoms  and  is 
expressed  as 

eT  *  Cv  (T  -  V  +  V  (8.2) 

To  determine  the  thermal  pressure  from  the 
thermal  energy,  eT.  the  first  law  of 

thermodynamics  and  a  Maxwell  relation  are  used  to 
oDtain 


«>T-  T  (ff)v-  P- 


(8.3) 

This  Is  a  general  equation  applicable  to  both 
components,  elastic  and  thermal.  Since  the  cold 
pressure  Is  only  a  function  of  volume  (or  more 
specifically  Interatomic  distance),  then  ?r  • 
Pc(v)  and  Eq.  (B . 3 )  gives  c 

de 

p  9  _ c 

c  “  ST 

aS  g1ven  Pr1or  t0  ^1s.  Eg. 

therma  I  'component!  a"aIyS,S-  F°r  the 

PT  ?  CvT  (^v  •  (8-4) 

Can  **  exPressed  In  terms  of  the 
coeff1*cfMt)Caser1Vat,Ve  r(v)  (the  Grune,sen 


C  T 

PT  T(v)-J- 


(8.5) 


where,  as  stated  earlier,  r(v)  «  y(|£.)  .  The 

measi!rpe^fC?hff,£1ent’  aS  exPr<?ss*P  here,  Is  a 
measure  of  the  change  In  thermal  pressure  to 
tnermai  energy. 

in  add???!!0?*  iMl*  (8*3)  40 P  (8.5)  are  solved 
In  addition  to  the  algebraic  component  relations. 


e  ■  ec  ♦  eT 

P  ■  pc  *  PT 

and  the  Hugonlot  Jump  conditions, 
1 


e  ■  e, 


!H  ■  2  VV  v) 


(8.6) 

(8.7) 

(8.8) 


Hgure*!?!*  ***  C°1d  c“"Press,°n  curve  shown  In 
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Fig.  11  Hugonlots  for  unreacted  (  *-0)  and  fully 
reacted  (  A *1 )  explosives 


For  this  particular  case,  the  ratio  r  /v  Is 
assumed  to  be  constant.  Figure  11  shows 
graphically  how  each  component  of  energy  Is 
represented.  The  area  A8D0A  Is  the  total 
increase  In  Internal  energy  due  to  shock 
pressure,  Pg.  Of  the  energy  deposited  by  the 
shock,  the  area  under  the  curve  FED  goes  into  the 
potential  or  elastic  energy.  The  remaining  area 
ABEFOA  Is  the  contribution  to  the  thermal 
component  and  therefore  determines  the 
temperature  rise  associated  with  the  shock. 
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